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We analyze greying of the dark soliton in a Bose- 
Einstein condensate in the hniit of weak interaction be- 
tween atoms. The condensate initially prepared in the 
excited dark soliton state is loosing atoms because of 
spontaneous quantum depletion. These atoms are de- 
pleted from the soliton state into single particle states 
with nonzero density in the notch of the soliton. As a 
result the image of the soliton is losing contrast. This 
quantum depletion mechanism is efficient even at zero 
temperature when a thermal cloud is absent. 



I. INTRODUCTION 

Since the experimental creation of atomic Bose- 
Einstein condensates (BEC) BEC is a subject of 
intensive experimental and theoretical research. The 
Gross-Pitaevskii equation (GPE) [|| is commonly used to 
describe with remarkable precision the ground state prop- 
erties and the dynamics of BEC. The GPE assumes that 
all atoms are in the same single particle wave-function. 
In this sense it is a time-dependent mean field equation. 
This equation is not an accurate description of the ground 
state in optical lattices. The experiment in the optical 
lattice trap ^ demonstrates substantial squeezing of the 
atomic ground state. As predicted in Ref. [Q, squeezing 
of the atomic state can also be achieved with a continu- 
ous phase contrast imaging of a condensate. The imaging 
is a non-destructive weak measurement of the density of 
atoms which makes the phase of the condensate wave- 
function ill defined. In this paper we ask a simple ques- 
tion: is the GPE an accurate description of a collectively 
excited condensate in a standard harmonic trap? 

The recent experiment by the Hannover group |^] 
strongly suggests a negative answer. In Ref. |^ exper- 
imental results are carefully compared with simulations 
of the GPE. The observed greying of the dark soliton is 
not consistent with the GPE. In this paper we analyze 
the problem in a weak particle interaction limit. The 
analysis allows us to relate greying of the dark soliton to 
a quantum depletion of atoms from the soliton state into 
single particle states with nonzero density in the soli- 
ton notch. This mechanism does not require any help 
from external agents like e.g. collisions with atoms from 
the thermal cloud |^ . The quantum depletion is present 
even at zero temperature. The condensate prepared in 
the excited dark soliton state has an intrinsic instability. 



Atoms coherently evolve away from the soliton conden- 
sate without help from any external dissipation. 



II. GROSS-PITAEVSKII EQUATION (GPE) 

N bosonic atoms with repulsive s-wave two-body in- 
teractions in a trap potential V{x) are described by the 
iV-body Schrodinger equation 
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* + a^(5(xfc -x,) ^ . (1) 



Here xi, . . . , xn) is a symmetric A^-body wave- 
function. We use a dimensionless form of the Schrodinger 
equation. When A'^ atoms are condensed in a single par- 
ticle state $(t,x) i.e. 



N 



*(i,xi,...,xjv) = [|$(i,xfc) , 



(2) 
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then evolution of $(t,x) can be described by the GPE 



-iv2 + y(x)+,g$*$ 



(3) 



with g = {N — l)a. The derivation of Eq.(||) from 
Eqs.(^||) assumes that the evolving state of A^ atoms 
remains in the subspace of product states (^). In this 
paper we show that this assumption is not satisfied for 
excited states of the GPE. Product states (^) are not 
eigenstates of exchange interactions which are neglected 
in the derivation of the GPE. 



A. Symmetries of the GPE 

GPE has symmetries similar to those of the one-body 
Schrodinger equation. For example, in one dimension an 
(anti-) symmetric ^{t,x), ^{t,—x) = {—)^{t,x), remains 
(anti-)symmetric in the course of the time evolution. In 
particular, a static dark soliton ||^ 



<I>(x) ~ tanh(Q;x) 



(4) 



in a symmetric trap remains antisymmetric and the den- 
sity of atoms n(x) — N\^{x)\'^ has a zero at x — Q. An 
antisymmetric dark soliton preserves its empty notch ac- 
cording to GPE. 
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III. THREE MODE APPROXIMATION FOR 
SMALL G 

We assume an effectively one dimensional cigar-shaped 
harmonic trap i.e. one with uj± ^ u>x when the transver- 
sal motion of atoms is frozen in the transversal ground 
state. The remaining x-dimension has the harmonic po- 
tential V{x) — Lo'^x^ /2. For small g we may truncate the 
single particle Hilbert space to the three lowest energy 
states (we choose oscillator units where lOx — ^) 



uo{x) 
ui{x) — X 

U2ix) 



^1/4 ' 



TT 



1/4 



(5) 



1*0 and 1*2 are symmetric and ui is antisymmetric. The 
(symmetric) ground state of the GPE is 



Mx) = Auoix) + (l-A2)i/2y2(x) , 



(6) 



with A E [0,1]. j4 = 1 for g = and decreases with 
increasing g. The admixture of U2{x) for nonzero g is the 
closest that we can get to the Thomas-Fermi limit in our 
3- mode approximation valid for small g. 

The antisymmetric excited state of the GPE 



^'soi(a;) = ui{x) 



(7) 



is a (small g) 3-mode analog of the dark soliton in Eq.(||). 

To study quantum evolution we will use the second 
quantized Hamiltonian 



H = Jdx ^dx^^dx^ + V{x)4>^4>+^^'i^'i^^ 



(8) 



In the 3-mode approximation ^(x) — aoUo{x)+aiUi{x) + 
02^2(2;), where Ui are usual bosonic annihilation opera- 
tors, and the Hamiltonian becomes 
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^ooooQ'o'^oQ'Oao + -^iiiiaiOiaiai + /2222a2'^2^2a2 
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-^0022 
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4aJd|aoai + ^ajajaidi + h.c.^ 
4dJdJdod2 + (^ajdja2d2 + h.c. 
4:a\a\aid2 + ^ajdjdidi + h.c. 
4djd|did2 + 2djd2didi + h.c. 
2d(')dJdod2 + h-C. 



+ -^^0222 200026202 -I- h.c 



(9) 



with integrals defined by 

/' + OC 
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dx Ua{x)ub{x)Uc{x)ud{x) 



(10) 



IV. DEPLETION OF THE GROUND STATE 

For noninteracting (g = 0) particle case modes (|^) con- 
stitute exact eigenstates of the system. If o^; 7^ there 
are various couplings between modes. The truncation 
of the Hilbert space to 3-modes only is a valid approxi- 
mation if the interaction energy per atom does not shift 
significantly energy of an atom, i.e. [N — l)/oooo < ^x- 
That allows us to employ values of g = {N — l)ax up to 
10. 

We have diagonalized for < 50 and for < 
[N — 1)03; < 5 in the Fock basis 



|7Vo,iVi,7V2) = 



(ot)^o(ot)^^(6t)^^ 
^/TVoWO/Val 



|0) 



(11) 



with No + Ni + N2^ N. We find that, except for [N - 
l)ax = 0, the iV-body ground state ICS*) is not exactly 
a product state like in Eqs.(^^. The second quantized 
version of that product state is 



\A) = 



^oj-f (l-^2)l/2^t 
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■10). (12) 



We find the product state which is the closest to the 
ground state \GS) with the help of the number of 
atoms in a given product state \A) defined by Na = 
{GS\d\dA\GS). The optimal product state \A) mini- 
mizes the fraction of atoms depleted from the condensate 



Dqs Min I 1 - 



Na 
N 



(13) 



The dependence of Dqs on the interaction strength g — 
{N — l)ax is shown in Fig.l. We find that A decreases 
with increasing [N — l)oa; (the condensate is transferred 
to the mode U2) and I?gs increases with [N — \)ax (in- 
creasing fraction of atoms is depleted from the conden- 
sate). As expected, see e.g. Ref. for the ground state 
of the system, the depletion is a small effect. 



V. DEPLETION FROM THE FIRST EXCITED 
STATE 

According to GPE the dark soliton in our 3-mode ap- 
proximation is the state (|7|j^), or in the second quantized 
language 



|0,7V,0). 



(14) 
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The state |0, A'', 0) is a stationary state of the GPE but 
it is not an eigenstate of H^. The quantum Hamilto- 
nian is allowed to mix antisymmetric modes with sym- 
metric ones provided the symmetry with respect to the 
parity operation in the full iV-particle Hilbert space is 
preserved. Various terms in H^, like e.g. a\a\aiai^ 
02020.101, or ojojoioi, are mixing the antisymmetric 
mode 1 with symmetric modes and 2. In more phys- 
ical terms, the vertex oJdjOiOi describes a collision of 
two atoms in the antisymmetric soliton mode 1 and their 
transfer into two symmetric modes and 2. The term 
a\a\aiai is again a collision of two atoms in the mode 1 
but this time both atoms are transferred to the symmet- 
ric mode 2. A similar vertex ojajoioi transfers pairs of 
atoms from the mode 1 to the symmetric mode 0. 



A. Degenerate subspace 

Among these processes the most efficient is 

h.c, (15) 



0^020.101 



because such a transfer preserves single particle energy. 
For ax = 0, the solitonic state |0,7V, 0) belongs to a de- 
generate subspace spanned by vectors 



\i,N - 2i,i). 



(16) 



When Qx 7^ 0, the process ( |l5| ) couples these states 
among each other. Due to the degeneracy, an arbitrary 
small Qx is enough to mix |0, N, 0) with other states from 
the subspace. Even for Oa, ^ all eigenstates of H3 in 
the subspace have nonzero depletion — even a very week 
interaction between atoms induces a significant depletion 
of the soliton. For the system in the ground state there is 
no similar process connecting degenerate states and the 
depletion in the ground state is very small. 

We have calculated the depletion of the eigenstates, 
\4>), of H3 from the soliton wave- function 



Ns 

Ds = l 



(17) 



where Ns — {(t)\a[ai\(t)) , for N = 30. In Fig. 2 we show 
the depletion of the eigenstate which has the smallest 
Ds as a function of {N — 1)02,. When {N — l)ax — *■ 0, 
the depletion does not tend to zero. This is because of 
the degeneracy of the subspace the soliton state belongs 
to for Qx — 0. Indeed, the |0,7V, 0) state is not a limit 
of any eigenstates of H3 when ax decreases. Such an 
effect is obviously possible for a harmonic potential only. 
For other traps, strong mixing of the soliton with the 
\i,N — 2i, i) states begins at finite o^,. The insert of Fig. 2 
compares the density profiles of the GPE state |0,30, 0) 
and of the least depleted 7?3-eigenstate. The eigenstate 
has a substantial nonzero density at a; = 0. 



B. Evolution from zero depletion state 

In contrast to the GPE, in the full quantum evolution 
there is no symmetry constraint to keep the notch of the 
initial soliton empty. As |0,iV, 0) is not an eigenstate of 
H^, the state of the system \il^{t)) prepared in the dark 
soliton state, |'!/'(0)) = |0, A^, 0), evolves away from this 
state. As a result of depletion the notch of the dark 
soliton is filling up and the soliton is greying. In Fig. 3 
we show the contrast C{t) as a function of time. The 
contrast is defined by 



C{t) = 



,{t) - n{t,x = 0) 
,{t) + n{t,x = 0)' 



(18) 



where nmax(i) is the maximal value of the density n{t, x) 
at time t. The insert of Fig. 3 shows the density n{x) of 
the state (initially prepared as the dark soliton) at dif- 
ferent moments of time. Due to the depletion the notch 
of the dark soliton is filling up. Imporantly, depletion of 
the soliton reaches much greater value than the deple- 
tion of the corresponding least depleted eigenstate of the 
system. 

The time-scale on which the soliton contrast decays 
can be estimated from below by A^/qooo- the interaction 
energy of atoms in the ground mode uq. /qooo > labcd so 
the interaction energy per atom in the A^-atom ground 
state can be used as an upper estimate for the deple- 
tion time from the soliton condensate. Employing this 
crude approximation to the experimental situation |^ 
(i.e. N = 1.5 ■ 10^ atoms of ^^Rb in the trap with 
iL;±_ — 271 X 425 Hz and lOx = 27r x 14 Hz) we get for 
the minimal depletion time tmin = 2Tr/D, = 0.6 ms. The 
interaction energy per atom, f2 = 11.2 kHz, has been es- 
timated employing the Thomas-Fermi approximation for 
a condensate wave-function. The minimal time imin, is 
about 15 times shorter than the observed experimentally 
greying time, i.e. ~ 10 ms |5j. We emphasize that this 
result is meant to give a scale on which the evolution of 
a soliton is affected by quantum depletion rather than a 
quantitative prediction for the experiment. More accu- 
rate estimate for the Thomas-Fermi limit of large g would 
have to include a large number of even and odd modes. 
In particular it would have to take into account depop- 
ulation of even modes by two-particle scattering to odd 
modes. Such calculations could in principle be done in 
the spirit of the formalism of Ref. |l^. However, in Ref. 
|lT| two of us present large g calculations in the frame- 
work of the Bogoliubov theory. They predict that a dark 
soliton will loose contrast after around 10 ms. 



VI. CONCLUSION 

We have analyzed quantum depletion of the solitonic 
solution of the GPE for weak interaction between atoms 
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[small g, see Eq. (^)]. We have shown that a condensate 
prepared initially in the soliton state loses its contrast 
according to the quantum multiparticle time evolution. 
Our 3-mode model is a suggestive example that the GPE 
cannot be used to describe large collective excitations of 
the condensate. The model is a proper description of the 
problem for a weak interaction between atoms only but it 
ilustrates essential processes responsible for the depiction 
of the solitonic wave-function. 

It was suggested in Ref. |^,|| that the deviations from 
GPE can be possibly explained by dissipation due to col- 
lisions of the condensate with the thermal cloud of non- 
condensed atoms. There is no doubt that dissipation 
influences dynamics of the soliton in that experiment. 
However, our analysis implies that a dark soliton may 
grey even in the absence of any thermal cloud. The notch 
may fill up with atoms quantum depleted from a conden- 
sate. 

Our analysis of the greying of the dark soliton may con- 
vey a wrong impression that a vortex in BEG should 
also grey as a result of quantum depletion. Solitons are 
experimentally excited by means of the so-called phase 
imprinting method In our (small g) 3-mode approx- 
imation we model the result of such an excitation as the 
product state of the antisymmetric modes Eqs. (0,^. In 
the case of a vortex state the situation is different. In- 
deed, vortices are created in a rotating (slightly asym- 
metric) cylindrical trap after achieving BEG or by cool- 
ing down the gas below the condensation temperature 
in an already rotating trap. Switching to the frame ro- 
tating with the trap one finds a single quantized vortex 
as the iV-atom ground state that is an eigenstate of the 
A^-atom Hamiltonian and as such it does not evolve in 
time. Eigenstates of a rotating condensate that consists 
of small number of atoms has been analyzed recently 
in Ref. It was observed that indeed the ground 

state with a vortex reveals a noticeable depletion. This 
ground state depletion is stationary, it does not grow in 
time. However, vortices created by the phase imprinting 
method or subject to suddenly stopped rotation are ex- 
pected to reveal much larger quantum depletion which 
may appear as greying. 

Note added: Quantum depletion in a stationary soli- 
tonic state has been predicted recently in . 
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FIG. 1. The fraction of atoms Dqs depleted from the con- 
densate in the A''-body ground state as a function of the inter- 
action strength (A'^ — l)ax for several values of A'' = 10, 30, 50 
(from top to bottom). 



FIG. 2. The depletion from the soliton wave-function 
[0,30,0) of the least depleted eigenstate of as a function 
of (N — l)ax- The insert shows the density of atoms n{x) 
in the dark soliton [0, 30, 0) (solid line) and in the eigenstate 
(dashed line) at {N - l)a^ = 4.7. 



FIG. 3. The contrast C{t) of the evolving soliton for 
(A'^ — l)ax = 1 and N = 20. The contrast will go through 
revivals for times later than 40 due to finite Hilbert space of 
the considered three mode approximation. The insert shows 
densities n{x) of the soliton at t = 10 (solid line) and at t = 30 
(dashed line). 
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